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Abstract. Based on detailed 2D numerical radiation hydro- 
dynamics (RHD) calculations of time-dependent compressible 
convection, we have studied the dynamics and thermal struc- 
ture of the convective surface layers of solar- type stars. The 
RHD models provide information about the convective effi- 
ciency in the superadiabatic region at the top of convective en- 
velopes and predict the asymptotic value of the entropy of the 
deep, adiabatically stratified layers (Fig. ||). This information is 
translated into an effective mixing-length parameter «mlt suit- 
able to construct standard stellar structure models. We validate 
the approach by a detailed comparison to helioseismic data. 

The grid of RHD models for solar metallicity comprises 
58 simulation runs with a helium abundance of Y = 0.28 in 
the range of effective temperatures 4300 K < T cS < 7100 K 
and gravities 2.54 < logg < 4.74. We find a moderate, never- 
theless significant variation of «mlt between about 1.3 for F- 
dwarfs and 1 .75 for K-subgiants with a dominant dependence 
on T e ff (Fig. ^|). In the close neighbourhood of the Sun we find 
a plateau where q;mlt remains almost constant. The internal 
accuracy of the calibration of cvmlt is estimated to be ±0.05 
with a possible systematic bias towards lower values. An anal- 
ogous calibration of the convection theory of Canuto & Mazz- 
itelli (1991, 1992; CMT) gives a different temperature depen- 
dence but a similar variation of the free parameter (Fig. ^). 

For the first time, values for the gravity-darkening expo- 
nent j3 are derived independently of mixing-length theory: 
(5 = 0.07... 0.10. 

We show that our findings are consistent with constraints 
from stellar stability considerations and provide compact fitting 
formulae for the calibrations. 
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1. Introduction 

The structure and evolution of late-type stars is intimately re- 
lated to convective transport processes taking place in their 
outer envelopes. Although the underlying physical principles 
are well known, the non-linear and non-local character of the 
equations describing the convective motions in a radiating, par- 
tially ionized fluid has hampered the development of a closed 
analytical theory. Hitherto standard stellar structure models 
include the convective energy transport in the framework of 
mixing-length theory (MLT). Despite its heuristic nature MLT 
has proven to be rather successful and is still the working horse 
in stellar structure modeling. During recent years much effort 
has been put into the refinement of the theoretical description 
of microphysical stellar plasma properties — the opacity and 
the equation of state. To fully benefit from these improvements, 
an accompanying development of the description of hydrody- 
namical transport properties appears to be necessary, in the first 
place aiming at a better understanding of the convective energy 
transport. In this paper we report on such an effort relying on 
direct numerical simulations of convective flows in solar-type 
stars. 

Primarily due to the heavy demands on computational re- 
sources, RHD simulations have concentrated very much on the 
Sun, and only rather few models have been constructed for 
other objects (Nordlund 1982, Steffen et al. 1989, Nordlund & 
Dravins 1990, Ludwig et al. 1994, Freytag et al. 1996). This 
sparseness has prevented a broader application of results of 
RHD simulations to the modeling of stellar atmospheres and 
evolution. We try to overcome this limitation in this work by 
computing a grid of RHD models, allowing a systematic study 
of the characteristics of convective flows in a larger region of 
the Hertzsprung-Russell diagram (HRD). We point out that the 
RHD simulations have now reached a level of sophistication far 
beyond idealized, simple numerical experiments. We consider 
it an important aspect of our work that it presents a quantita- 
tive study and — as such — puts emphasis on the estimation of 
the achieved accuracy. Aiming at an application in the field of 
stellar evolution, we concentrate on the evaluation of one struc- 
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tural quantity, namely the entropy of the adiabatically stratified 
layers deep in the convective stellar envelope. In the following 
we shall refer to this entropy as s cnv . Its relevance for stellar 
structure stems from the fact that it strongly influences the ra- 
dius of a star. To obtain handy numbers, all entropy values in 
this paper — unless stated otherwise — are given in units of 
So = 10 9 erg/g/K. To link our results more directly to standard 
stellar structure modeling we translate s cnv into an equivalent 
mixing-length parameter. 

Some words concerning our nomenclature: we use the 
term "solar-type" for stars with extended convective envelopes 
where the thickness of the superadiabatic layers at the top 
of this envelope is small in comparison to the stellar radius. 
"Grey" radiative transfer means that frequency-independent 
(mean) opacities were used in the computation of the radiation 
field. The opacities still depend on temperature and density and 
include contributions from spectral lines. 

In the paper we proceed as follows: we start with method- 
ical aspects describing our hydrodynamical models, the basic 
idea and the procedure to derive s cnv , and the translation of s en v 
into an equivalent mixing-length. We continue with the valida- 
tion of our method by showing that we are able to predict the 
solar structure derived from helioseismic measurements within 
small uncertainties. We then present the calibrations of MLT 
and CMT We discuss the application of our results to stellar 
modeling, point out consequences of stellar stability consider- 
ations, present a derivation of the gravity-darkening exponent, 
and contrast our approach with others. We conclude with fu- 
ture perspectives. In the appendix we provide some auxiliary 
data helping to utilize our findings in stellar structure models. 

2. Methodical aspects 

2.1. Hydrodynamical models of solar-type surface convection 

We have obtained detailed 2-dimensional models of the sur- 
face layers of solar-type stars from extensive numerical simu- 
lations solving the time-dependent, non-linear equations of hy- 
drodynamics for a stratified compressible fluid. The calcula- 
tions take into account a realistic equation-of-state (EOS, in- 
cluding the ionization of H and He as well as formation of 
H2 -molecules) and use an elaborate scheme to describe multi- 
dimensional, non-local, frequency-dependent radiative trans- 
fer. Similar to classical model atmospheres, the hydrodynami- 
cal models are characterized by effective temperature T c g, ac- 
celeration of gravity log g, and chemical composition. They in- 
clude the photosphere as well as part of the subphotospheric 
layers, with an open lower boundary, allowing a free flow of 
gas out of and into the model. A fixed specific entropy s* is 
(asymptotically) assigned to the gas entering the simulation 
volume from below. The value adopted for s* uniquely deter- 
mines the effective temperature of the hydrodynamical model. 
For details about the physical assumptions, numerical method 
and characteristics of the resulting convective flows see Ludwig 
et al. (1994) and Freytag et al. (1996). 



2.2. From the surface to the base of a convection zone 

Figure [j]a shows the mean entropy as a function of depth ob- 
tained from a hydrodynamical granulation model of the Sun 
by averaging over horizontal planes and over time. As in this 
example, our models in general do not extend deep enough to 
include those layers where the mean stratification of the con- 
vection zone becomes adiabatic. While the mean entropy strat- 
ification of the hydrodynamical models does not permit a direct 
determination of the entropy corresponding to the adiabat of the 
deep convection zone, the spatially resolved entropy profiles 
contain additional information. Figure |l]b displays the entropy 
profiles for an arbitrary instant of the sequence from which the 
mean stratification in Fig. [j]a was computed. The granular con- 
vection pattern at the surface of solar-type stars is formed by 
broad hot upflows accompanied by concentrated cool down- 
drafts. Figure [j]b shows a remarkable entropy plateau in the 
subsurface layers, indicating that — in contrast to the narrow 
downdrafts — the gas in the central regions of the broad as- 
cending flows is still thermally isolated from its surroundings. 
Neither radiative losses nor entrainment by material of low en- 
tropy can produce significant deviations from adiabatic expan- 
sion until immediately below the radiating surface layers. The 
height of the entropy plateau is essentially independent of time 
and corresponds to s*. 

We suggest that s* may be identified directly with the en- 
tropy Senv °f the deep, adiabatic convective layers: 

Scnv — S . ( 1 ) 

This idea has been put forward by Steffen (1993) and Ludwig 
et al. (1997) and is based on the qualitative picture of solar- 
type convection zones proposed by Stein & Nordlund (1989) 
(hereafter referred to as "SN scenario") which is fundamen- 
tally different from MLT assumptions. According to this sce- 
nario the downdrafts continue all the way from the surface 
to the bottom of the convection zone, merging into fewer and 
stronger currents at successively deeper levels. The flow closes 
only near the base of the convective envelope. Most of the gas 
elements starting from the bottom of a deep convection zone 
overturn into neighbouring downflows before reaching the sur- 
face. Only a very small fraction of gas continues to the surface, 
reaching the layers corresponding to the location of the lower 
boundary of our hydrodynamical models essentially without 
entropy losses, following an adiabat almost up to the visible 
surface. Hence, s* obtained from the simulations is the entropy 
of the warm, ascending gas throughout the convection zone. 
This, in turn, is very nearly equal to the mean (horizontally 
averaged) entropy s cnv near the base of the convection zone 
because (i) the downflows are markedly entropy-deficient only 
near the surface and become continuously diluted by overturn- 
ing entropy-neutral gas as they reach greater depths, and (ii) 
the fractional area occupied by the downdrafts decreases with 
depth. 

We note that for this investigation it is actually irrelevant 
whether the downflows possess a plume-like character or take 
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Fig. 1. Depth dependence of the entropy in the solar surface layers as obtained from hydrodynamical simulations (T c g = 5770 K, 
log g = 4.44, model code L71D07) performed on a 140 (x) by 71 (z) grid with frequency-dependent radiative transfer. The mean 
entropy (horizontal and temporal average) is shown in panel (a), spatially resolved entropy profiles in panel (b). Geometrical 
height zero corresponds to tr ss = 1 ■ Note that the model comprises only the uppermost part of the 200 Mm deep solar convective 
zone. 



place in another form. The essential point is that sinking ma- 
terial does not locally affect the entropy of rising material by 
draining heat from it. The role of the downflows is reduced to a 
dynamical one: while sinking downwards they simply displace 
buoyancy-neutral material and push it upwards. 

2.3. Functional dependence of s onv on stellar parameters 

When considering the MLT picture and the SN scenario it is 
interesting to ask what stellar properties determine s cnv or — 
in other words — what coordinates are appropriate to describe 
its functional dependence in the HRD. We argued that hydro- 
dynamical models of the surface layers are able to provide in- 
formation about Scnv The models are characterized by the at- 
mospheric parameters and consequently we describe s cnv in 
terms of them. Whether the standard atmospheric parameters 
T e g and \ogg are the most suitable coordinates is not clear. 
One might speculate that e.g. the surface opacity is a physi- 
cally more relevant quantity. Nevertheless, for solar-type stars 
the conditions at the stellar surface govern the global envelope 
structure and the standard atmospheric parameters are suitable 
coordinates to parameterize them. They have the advantage that 
they are external control parameters and not part of the solution 
of the problem. For solar-type stars with superadiabatic regions 
which are thin in comparison to the stellar radius we expect the 
same qualitative behaviour irrespective of whether we consider 
the MLT picture or the SN scenario. Global stellar parameters 
(mass, radius, or age) play only an indirect role for the entropy 
jump. The situation changes when the size of the granular cells 
or the thickness of the superadiabatic layer become comparable 
to the stellar radius. This might happen in giants and one has 
to account for effects introduced by the global stellar structure, 
i.e. sphericity effects. 



2.4. ckmlt from envelope models 

Although Scnv allows the construction of the envelope structure 
with the necessary precision, it is of interest to translate this 
quantity to an equivalent auur since this parameter is con- 
ventionally used in stellar structure models. For this purpose 
we computed grids of standard stellar envelope models (not 
subject to central boundary conditions) based on MLT, cov- 
ering the relevant range of effective temperature and gravity 
and assuming Y = 0.28 and solar metallicity for the chemi- 
cal composition. In these models «mlt was a free parameter. 
For given stellar parameters (T c g , log g) we obtained s cnv as 
a function of aiviLT- By matching s on v from the RHD models 
we deduced the corresponding «mlt- Figure ^ illustrates this 
procedure for a number of representative models. Each panel 
shows the entropy profile of the envelope model matching s onv 
from a RHD simulation. For comparison two further envelope 
models are plotted with «mlt varied by ±0.05 with respect to 
the matching one, as well as an estimate of the uncertainty of 
Scnv stemming from the temporal fluctuations of s* which we 
observe during the simulation run. Although at first glance the 
matching is a straightforward procedure, it is important to take 
care of a number of fine points in order to ensure a unique and 
well-defined calibration of q;mlt- We discuss these points in 
the following. 

Concerning the opacities and EOS we tried to stay as close 
as possible to the physical description used in the RHD models. 
Due to this differential procedure, systematic uncertainties are 
substantially reduced. In the RHD models ATLAS6 opacities 
according Kurucz (1979) are adopted. These opacities contain 
no contribution from molecular lines. Since the envelope mod- 
els reach much deeper than the RHD models we supplemented 
the ATLAS6 opacities by OPAL opacities for higher tempera- 
tures (Rogers & Iglesias 1992). Care was taken that the chem- 
ical mixture assumed in both data sets was as similar as possi- 
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ble. In the envelope models we used the RHD EOS throughout. 
It is a simple EOS taking into account the ionization equilib- 
ria of hydrogen and helium according to the Saha-Boltzmann 
equations. 

The RHD EOS has the advantage of great smoothness and 
computational efficiency. It might appear that this description is 
too crude, especially for the envelopes of cooler objects where 
non-ideal effects become more pronounced. However, this sim- 
plification is tolerable since the entropy profile, which is the 
major concern here, is not sensitively dependent on the EOS. 
The entropy gradient is proportional to the difference between 
actual and adiabatic temperature gradient which in MET is re- 
lated to the convective flux. Hence, the constraint that a nomi- 
nal total flux has to be transported essentially fixes the entropy 
gradient irrespective of the EOS employed, while the tempera- 
ture profile is noticeably affected by the choice of the EOS. 

Following the conventional way to treat the atmosphere in 
stellar structure models we prescribe a certain T(r)-relation 
representing the temperature run in the optically thin regions. 
Its significance for the envelope structure emerges from the 
fact that it determines the level of the entropy minimum in the 
deeper atmosphere where the stratification becomes convec- 
tively unstable (the starting level for the entropy jump). In evo- 
lutionary calculations various descriptions of the atmosphere 
are commonly used, differing in the level of sophistication be- 
tween simple analytical or empirical T(r)-relations and full- 
fledged model atmospheres. Here we chose a T(r)-relation that 
mimics closely the average atmospheric structure of the hydro- 
dynamical models. We reproduce in particular the atmospheric 
entropy minimum found in the hydrodynamical models. The 
hydrodynamical models predict an average temperature struc- 
ture in the deeper photosphere (around — 2 < logr ross < 0) 
that closely resembles a stratification in radiative equilibrium. 
We emphasize that this statement refers to an averaging of the 
temperature on surfaces of constant optical depth in the RHD 
models. This averaging procedure preserves the radiative flux 
properties of the RHD models to a good approximation (see 
Steffen et al. 1995), and hence appears suitable to characterize 
their atmospheric structure in the present context since radia- 
tive surface cooling plays the dominant role for the convection 
driving. 

In practice, we used an analytical fit to the exact 
grey T(r)-relation for an atmosphere in radiative equilib- 
rium according Unsold (1955) in envelope models intended 
for comparison with RHD models adopting grey radia- 
tive transfer. From the small set of RHD models adopting 
frequency-dependent radiative transfer we derived the average 
T(r)-relation and constructed an analytical fit to it. Later in 
this paper we shall discuss the role of grey versus frequency- 
dependent radiative transport in the RHD models, demonstrat- 
ing the need for a clean distinction between both types of 
T(r)-relations. 

The influence of the T(r)-relation grows as the entropy 
jump itself decreases, so it is largest for models at low ef- 
fective temperature and high gravitational acceleration. In ear- 
lier papers on the subject (Ludwig et al. 1996, 1997) we used 



T(r)-relations which were not selected on a strictly differen- 
tial basis, leading to systematic differences between the older 
and the present calibration of «mlt- Figure [| illustrates the 
match of the atmospheric entropy minima. (The range in opti- 
cal depth over which the RHD models are plotted do not rep- 
resent the full extension of the RHD models.) Our present pro- 
cedure gives now an improved match with the largest deviation 
found for the (4500, 4.44) model. We expect that a further im- 
proved fit of the atmospheric stratification of the RHD model 
would lead to an increase of Q!mlt- But note that using the 
ctMLT derived from models in this region of the HRD would 
produce a consistent s onv in a stellar structure model as long 
as it is basing on a description of the stellar atmosphere by 
a T(r)-relation resembling an atmosphere in radiative equilib- 
rium. 

As a side-issue the Fig. ^demonstrates that in RHD models 
the higher atmospheric layers are systematically cooler than in 
an atmosphere in radiative equilibrium. This is due to the sit- 
uation that in a convective atmosphere the temperature is gov- 
erned by the competition between adiabatic cooling of over- 
shooting material and heating by the radiation field. Quan- 
tatively, our grey RHD models overestimate the temperature 
reduction since they do not include the important heating by 
spectral lines in the higher layers. However, as we will see in 
the next section this is of minor importance for the determina- 
tion of aMLT- The deviations between average RHD stratifica- 
tion and envelope models underline that our calibrated qmlt 
only fits the asymptotic entropy and the entropy jump. The 
overall stratification is not well matched by an envelope model 
with the fitted q;mlt- 

Effects of turbulent pressure are ignored in the envelope 
models. This is not only done for computational convenience 
but also reflects our opinion that the recipes provided for its 
treatment in the framework of MET overestimate the effects, 
doing more harm then good when included. In MLT the con- 
vective motions are restricted to the unstable layers according 
to the Schwarzschild criterion. At the upper boundary of a sur- 
face convective zone this produces a sharp decline of the con- 
vective velocity amplitude and a correspondingly large gradient 
of the turbulent pressure that alters the hydrostatic structure of 
a model significantly. In contrast, hydrodynamical models (see 
Freytag et al. 1996) predict a rather smooth run of the velocity 
amplitude in that region with a correspondingly smaller gradi- 
ent of the turbulent pressure. Clearly, the inclusion of turbulent 
pressure in ID envelope models requires a proper treatment of 
overshooting. 

Finally, it is important to realize that a unique formulation 
of MLT does not exist. Rather, different versions of MLT are in 
use. In the following we refer to the version of MLT originally 
given by Bohm-Vitense (1958). To further elaborate this point, 
we explicitly present details of our MLT implementation to- 
gether with commonly used other formulations in Appendix A. 
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Fig. 2. Representative examples of entropy profiles (in units of 
s = 10 9 erg/g/K) as a function of Rosseland optical depth 
logr ross . Each panel shows the r-averaged entropy of a RHD 
model (thick solid line) with envelope models (thin solid lines) 
for the matching Qmlt and «mlt ± 0.05 for a given T c g/ log g 
combination. In all models the radiative transfer was treated in 
the grey approximation. The dotted horizontal lines indicate s* 
of the RHD model and its temporal RMS-fluctuations. Note the 
significantly different entropy jumps. 

3. The solar benchmark 

The natural benchmark for the scenario described above is of 
course the Sun. Not only are its global stellar parameters known 
with exceptional accuracy, but helioseismology has provided us 
with detailed information about the structure of the solar inte- 
rior. Recent astrophysical convection theories were developed 
in the context of stellar evolution and tests of these theories 
were devised within that framework, in particular by studying 
their effects on the evolution of the Sun (Lydon et al. 1992, 
1993a; Canuto & Mazzitelli 1991, 1992; Canuto et al. 1996). 
In our opinion, a more direct way for a validation of such the- 
ories is to look at the present solar structure without taking re- 
course to evolutionary calculations. This avoids the possibility 



Table 1. Results from hydrodynamical (RHD) models and 
helioseismology for different helium abundances. Y: helium 
abundance, s cnv : entropy in the adiabatically stratified layers in 
units of So = 10 9 erg/g/K, TV,,: number of frequency bands for 
modeling the radiative transfer, cxmlt'- mixing-length parame- 
ter. The data for Y = 0.246 are obtained by linear interpolation 
between the data for Y = 0.24 and 0.28. 



Source 


Y 


N v 


Scnv / S() 


QMLT 


RHD 


0.28 


1 


1.723 


1.577 


RHD 


0.28 


5 


1.738 


1.588 


RHD 


0.24 


1 


1.782 


1.609 


RHD (inferred) 


0.24 


5 


1.800 


1.609 


Helioseismology 


0.24 




1.778 


1.727 


Helioseismology 


0.28 




1.736 


1.594 


RHD 


0.246 


5 


1.791 ±0.01 


1.61 =F 0.05 


Helioseismology 


0.246 




1.771 ± 0.005 


1.71 =F 0.02 



that evolutionary changes of the solar structure interfere with 
model properties that are used for the assessment of the accu- 
racy of the convection theory under consideration. 

In the following we shall compare our RHD model predic- 
tions for s cnv with related helioseismic measurements. On an 
absolute scale any deviation from the actual physical situation, 
e.g. in composition, shows up as mismatch between observa- 
tions and theoretical predictions without being necessarily re- 
lated to flaws in the convection model. In our solar models we 
paid particular attention to the role of the helium content Y and 
details of the radiative transfer. 

Table 1 summarizes the findings for the Sun. For inter- 
preting the different entropies given in the table we note here 
that the entropy jump over the superadiabtic layers in the Sun 
amounts to « 0.2 So- The first block of data (4 entries) in Ta- 
ble 1 refers to RHD models. For solar effective temperature 
and gravitational acceleration we calculated RHD models with 
Y = 0.24 in addition to models with our standard Y = 0.28. 
We considered RHD models where the radiative transfer was 
treated in grey approximation (N v = 1), and where its fre- 
quency dependence was approximately taken into account us- 
ing 5 representative frequency bands {N v — 5, for details about 
the radiative transfer scheme see Ludwig et al. 1994). The en- 
tropies Senv of the adiabatic part of the convective envelopes are 
based on the assumption s cnv = s*. The entropies derived from 
the RHD simulations were converted into an effective mixing- 
length parameter cumlt as outlined above. 

There is a noticeable dependence of s onv on the treatment 
of the radiative transfer and the helium content. The different 
treatment of the radiative transfer leads to a different atmo- 
spheric temperature structure, resulting in a change of the ab- 
solute entropy level in the deeper atmosphere. There is some 
further influence on the entropy jump itself since the radiative 
energy exchange in the surface layers is slightly modified. But 
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despite significant differences in s cnv for the two RHD models 
with Y = 0.28, both models give the same cumlt values within 
the uncertainty of q:mlt of ±0.05. This is achieved by com- 
paring RHD models and envelope models on a strictly differ- 
ential basis. The T(Y)-relation in the envelope models was se- 
lected to closely follow the relation found in the RHD models: 
RHD models basing on grey radiative transfer were compared 
to envelope models calculated with a grey T(T)-relation in their 
atmospheric layers. For RHD models basing on frequency- 
dependent radiative transfer, a T(r)-relation was derived by 
horizontally and temporally averaging their temperature struc- 
ture on surfaces of constant optical depth. This T(T)-relation 
was subsequently used in the calculation of the related en- 
velope models. When mixing solar RHD and envelope mod- 
els with incompatible T(r) -relations, «mlt changes by up 
to 0.1. We have performed the same tests for a hotter RHD 
(T c ff = 6500 K) model with similar findings. The insensitivity 
of q;mlt with respect to the treatment of the radiative transfer 
allowed us to largely rely on computationally less demanding 
RHD models employing grey radiative transfer for studying the 
scaling of «mlt across the HRD. 

The entropy difference due to a change in the helium con- 
tent is primarily related to the associated change of the mean 
molecular weight. To first order, the RHD models react to 
changes in Y by keeping the temperature-pressure structure 
invariant. The density scales in proportion to the change in 
molecular weight (here: fi(Y = 0.28)//x(Y = 0.24) = 1.037); 
entropy and internal energy per unit volume remain the same, 
allowing the model to transport the same energy flux without a 
change of the flow velocity. The simple scaling behaviour of the 
RHD models with Y and the invariance of ckmlt found in the 
case Y = 0.28 lead us to derive the values for Y = 0.24, N„ = 
5 from an envelope model assuming the same «mlt as in the 
case Y = 0.24, N„ = 1. The entry for this model in the table is 
flagged as "inferred". 

The second block of data (2 entries) in Table 1 shows re- 
sults derived from helioseismic considerations. Helioseismol- 
ogy provides the sound speed-density profile in the Sun with 
high accuracy. In the deeper, adiabatically stratified layers of 
the envelope this profile defines an adiabat in the sound speed- 
density plane. If chemical composition and the EOS are known 
one can label this adiabat with its actual entropy value. Note 
that only thermodynamic properties enter into this procedure; it 
is independent of opacity or results from stellar structure mod- 
els as long as one regards the helioseismic measurements them- 
selves as independent. The helioseismic entropy values of the 
second block of data in Table 1 are derived in this way by as- 
suming a certain Y and adopting the OPAL EOS (Rogers et al. 
1996) to establish the relation between sound speed, density, 
and entropy. The «mlt values are again derived from enve- 
lope models that match s onv . In contrast to s onv , the determina- 
tion of q;mlt depends on low-temperature opacities and stellar 
structure considerations. The helioseismic data are due to Antia 
(1996) and were kindly made available to us by the author. 

The third block of data (2 entries) finally summarizes our 
best estimates for s cnv from RHD models and helioseismology. 



They were obtained by linear interpolation between the data 
for Y = 0.24 and Y = 0.28 to Antia's (1996) best helioseis- 
mic estimate for the helium content Y = 0.246 based on the 
OPAL EOS. The uncertainties given in the helioseismic case 
are dominated by the uncertainty of the helium content. We 
have adopted an uncertainty of ±0.005, roughly comprising 
the range of Y discussed in the literature. In the RHD case the 
uncertainties reflect the variations which we find between sim- 
ulation runs for the same physical parameters that differ in nu- 
merical details (grid resolution, size of the computational box, 
explicit or implicit treatment of the energy equation, etc.). 

Of course, all our error estimates do not include possible 
systematic effects. The helioseismic Y can be affected by sys- 
tematic errors in the EOS. The RHD results are certainly af- 
fected by the two-dimensionality of the models. Indeed, we 
interpret the discrepancy between RHD and helioseismic re- 
sults as systematic effect: our 2D RHD models overestimate 
Sonv by 0.02 s (10% of the overall entropy jump between the 
surface and the adiabatic layers) and correspondingly underes- 
timate the solar o;mlt by 0.1. This conclusion is supported by 
first results from a differential comparison of 2D and 3D RHD 
simulations for the Sun. The 3D run predicts a smaller s en v by 
0.013s , corresponding to a larger cumlt by 0.07. Taking this 
fact into consideration, the theoretical and observational deter- 
mination of ttMLT become fully consistent. 

We expect another systematic effect to influence the abso- 
lute value of q;mlt affecting both the observationally and the 
theoretically derived values in the same manner. There exist in- 
dications that our low-temperature opacities underestimate the 
real stellar opacity in the deeper, subphotospheric layers. Our 
differential approach for determining cumlt cannot fully com- 
pensate for this deficiency since the opacities affect RHD and 
MET models differently. In RHD models the opacity effects 
are confined to the very thin cooling layer on top of upflowing 
regions, while in MET models the opacities are important all 
over the zone of high superadiabaticity. This biases our cumlt 
values again towards lower values. In the region of high supera- 
diabaticity our opacities are presumably too small by 10-20%, 
causing cumlt to be underestimated by « 0.1. 

However, in view of the fact that the RHD models are es- 
sentially parameter-free we consider the match in 2D as satis- 
factory. Besides the zero point, the functional dependence of 
Scnv and aMLT on the fundamental stellar parameters is of ma- 
jor interest. We believe that for these scaling relations our error 
estimates based on purely intrinsic model uncertainties are well 
justified. 

4. Calibration data for solar-type stars 

As for the Sun, we studied the behaviour of s cnv and ckmlt 
over the HRD in the solar neighbourhood. This investigation 
demanded a large number of simulation runs since we had to 
cover many objects and had to clarify how densely the HRD 
should be sampled in order to record the relevant variations of 
c*mlt- The results are now based on 58 simulation runs for 
solar metallicity and Y = 0.28, covering the range 4300 K < 
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T eS < 7100 K and 2.54 < log 5 < 4.74. Corresponding results 
for metallicities down to [M/H] = — 2 will be published in a 
subsequent paper. 

Various aspects determined the region of the HRD which 
is covered by our models: Towards higher T c s the convective 
zones become shallow and hence insignificant for the over- 
all stellar structure. Towards lower T c g molecular opacities 
become more and more important; these are not included in 
our opacities yet. Moreover, convection becomes very efficient, 
leading to small entropy jumps which cause problems in the nu- 
merical modeling. Towards lower log g the large entropy jumps 
cause other numerical problems related to the steep tempera- 
ture and velocity gradients in the surface cooling layer. And 
finally, no solar-type stars are located towards higher log g. 

We present our findings in graphical as well as numerical 
form. Like in the solar case, several RHD simulation runs dif- 
fering in numerical details were performed for certain atmo- 
spheric parameter combinations, thus allowing us to estimate 
the internal uncertainties associated with our models. In the 
figures we provide numerical values and show polynomial or 
exponential least-squares fits to the data points. The fits are 
constructed such that the residuals between fitted and actual 
data values show no remaining systematic dependence on the 
stellar parameters. In Appendix C we give explicit expressions 
for these fits. In this way our results can be easily utilized 
in other applications. We made efforts to produce fits with a 
simple and smooth functional dependence. Nevertheless we re- 
mind the reader to the fact that polynomial fits of higher order 
are not well suited for extrapolation. In other words: the fits 
quickly loose their meaning outside the region covered by the 
grid of RHD models and should not be extrapolated too far! 

Figure || displays s cnv as a function of T e g and log g in the 
HRD. In the region where s cnv > 2.6 x 10 9 erg/g/K, convection 
becomes inefficient and the convective zones are superadiabatic 
throughout. Entropies in that region refer to the value encoun- 
tered at the base of the convective envelope and are still derived 
assuming s cnv = s*. Besides s cnv itself, the entropy jumps are 
of interest since they measure the efficiency of the convective 
energy transport. The entropy jumps shown in Fig. |] are the en- 
tropy differences between the photospheric entropy minimum 
on the optical depth scale and s cnv (see Fig. ||). 

Figure [| shows that the RHD model grid covers more than 
an order of magnitude in the entropy jump, comprising stars 
with rather efficient envelope convection (with small entropy 
jumps) in the regions of K-dwarfs and stars with inefficient 
convection (large entropy jumps) in the region of F-dwarfs. 
The significant variation of the entropy jumps indicates that the 
simple dependence of s on v seen in Fig. g is non-trivial. It is not 
just reflecting the changes of the stellar surface conditions due 
to changes of opacity and ionization balance but carries infor- 
mation about changes in the efficiency of the convective energy 
transport. 

With the aid of envelope models (as described before) 
we translated the s cnv values given in Fig. |3| into equivalent 
mixing-length parameters which are shown in Fig. ||. We find a 
moderate, nevertheless significant variation of cumlt between 



about 1.3 for F-dwarfs and 1.75 for K-subgiants. In the close 
neighborhood of the Sun we find a plateau where q;mlt re- 
mains almost constant. As discussed previously, the absolute 
values of «mlt are probably less reliable than its scaling prop- 
erties as displayed in Fig. |[ 

The theory of Canuto & Mazzitelli (1991,1992; CMT) de- 
scribes convection analytically within the picture of a turbu- 
lent medium. In view of the significant amount of work that 
went into the implementation and verification of this prescrip- 
tion by several groups, it appeared worthwhile to compare our 
results with CMT. Like MLT, CMT contains a free length-scale 
A. According to Canuto & Mazzitelli A is essentially the dis- 
tance to the upper (Schwarzschild) boundary of the convective 
envelope; for purposes of fine tuning they add a small fraction 
of the pressure scale height at the upper boundary. We follow 
this recipe here by writing A = z + «cmt Hp t0 p and trans- 
lated our s cnv to acMT values, again with the help of enve- 
lope models, but now based on CMT. As shown in Fig. ^, we 
find acMT ~ 0.4 for the Sun. To assess the sensitivity of s en v 
to acMT, we note that the overall entropy jump at fixed solar 
T e g and log g increases by roughly a factor of 2 if one reduces 
acMT from 0.4 to zero. As in the MLT case, we find a region 
of rather constant acMT in vicinity of the Sun and a significant 
variation between zero in the region of K-dwarfs and almost 
0.6 for F-dwarfs. Interestingly, for the model with the lowest 
T e g, we find an «cmt slightly smaller than zero. This value 
was formally derived by extrapolating the s en v(c*CMT)-relation 
towards negative acMT- 

5. Implications and discussion 

5.1. Calibrating cumlt with evolutionary models of the Sun? 

The RHD models predict almost no change of cumlt during 
the solar main-sequence evolution and at most a very weak de- 
pendence of the convective efficiency on the helium abundance 
in the convective envelope (see Fig. || and Table 1). This is an 
important piece of information since it justifies the common 
practice of calibrating the solar cumlt with evolutionary mod- 
els. However, we think that the accurate knowledge about the 
solar envelope structure from helioseismology now provides a 
cleaner way to calibrate the efficiency of solar convection by 
the procedure outlined in Sect. ||. Calibrating q;mlt with the 
help of the present envelope structure would provide insight 
into the validity of parts of the physics put into the evolution- 
ary model calculations. 

For instance, calibrated solar models with and without ele- 
ment diffusion, respectively, provide different values for Qmlt 
(cf. e.g. Richard et al. 1996). Models including diffusion give 
a lower helium abundance at the surface for the present Sun 
and a somewhat higher q;mlt than models without diffusion. 
This might appear to contradict our findings, namely that «mlt 
is almost independent of the helium content. The reason is 
of course that the evolutionary models without diffusion do 
not represent the Sun since their present surface helium abun- 
dance is not the actual one. This shortcoming of an evolution- 
ary model without diffusion is compensated for by a change of 
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Fig. 3. s C nv in units of s = 10 9 erg/g/K from the grid of RHD models for Y = 0.28 and solar metallicity. Symbols indicate RHD 
models. Attached to the symbols the actual s onv values are given; the contour lines present a fit to them. The fitting function is 
given in Appendix C. 



Q;mlt- As a price, no assessment of the quality of the model 
or the adequacy of the input physics is possible any more. If, 
however, one calibrated cumlt at the present envelope struc- 
ture, shortcomings in the input physics would not be hidden 
in the auirr calibration but would show up as a mismatch be- 
tween the predicted and the actual solar radius. 

5.2. How to make best use of the calibration data in 
evolutionary models? 

In the case of the Sun our calibration underestimates the abso- 
lute value of the mixing-length parameter by about 0.1 .. . 0.2, 
a value which can be explained by a combined systematic ef- 
fect due to the too small low-temperature opacities and the 2D 
approximation. The scaling behaviour of the «mlt with ef- 
fective temperature and surface gravity is a differential result 
and will be much less affected by these systematic shortcom- 
ings. From this perspective it seems that a scaling of the «mlt 
values presented here by a constant factor (slightly larger than 
unity) is still admissible. In the context of stellar evolutionary 
models we suggest to calibrate omlt at the present Sun and 
use the ratio to qmlt from our calibration as a scaling factor. 
We propose a constant scaling factor since it appears plausi- 
ble to compensate for the systematic offset seen in the Sun in 
a simple, homogeneous fashion. For the time being this is the 



best one can do since the Sun is the only star where we know 
ttMLT with sufficient accuracy to critically test our results ob- 
servationally. However, we emphasize that at the moment the 
scaling is a well-motivated modification to our results. We ex- 
pect to eliminate the necessity of this step by improving on the 
identified systematic shortcomings of our approach in future 
models. 

The simplified radiative transfer in the RHD models pro- 
vides only an approximate atmospheric T(r)-relation. We have 
seen by comparison of grey and non-grey models that the in- 
fluence of the T(r)-relation on the determination of cumlt can 
be strongly reduced if one sticks to the same approximation 
in the envelope and RHD models. In this way we separated 
to some extend properties of convection from issues related to 
opacities and radiative transfer. If one now aims at a compari- 
son with observations, the radiative properties should be mod- 
eled as accurately as possible. In evolutionary models intended 
to reproduce observations (e.g. colors) one should use as ac- 
curate T(r)-relations as possible, preferentially derived from 
full-fledged model atmospheres. The qimlt values from our 
calibration are not affected by such a change to a moderately 
differing T(r)-relation. 



Ludwig et al.: A calibration of the mixing-length for solar- type stars I. 9 




7000 6500 6000 5500 5000 4500 

T eff [K] 



Fig. 4. Entropy jumps in units of Sq = 1 9 erg/g/K. The entropy jump was calculated as difference between s cnv and the entropy 
at the photospheric entropy minimum. See Fig. || for further explanations. 



5.3. Consistency with stellar stability 

A variation of «mlt and corresponding variation of s cnv with 
the atmospheric parameters T c ff and log g can lead to a secular 
instability of a star on the Kelvin-Helmholtz timescale of its 
convective envelope. In order to see this, let's consider a small 
perturbation of s cnv in a star's envelope. As known from stel- 
lar structure models such a disturbance alters its radius while 
the luminosity remains largely unchanged. The change in ra- 
dius leads to a change of the surface parameters. Following 
Christensen-Dalsgaard (1997, his relations (10) and (12)) one 
approximately finds for the differential relation between s cnv 
and T c ff valid for stars similar to the Sun 



<9Sei 

9 T ff J L=CO nst 



T 



cir 



(2) 



where c p is the typical specific heat at constant pressure in the 
convective envelope. 

While relation (Q) describes the response of the global stel- 
lar structure on changes of s cnv the question is whether this 
is compatible with the surface conditions. Since from the sur- 
face conditions s onv is a function of T e fj and log g, there is the 
possibility of a feedback on s onv which can amplify or damp 
an perturbation of s en v One would obtain an amplification if 
Sonv drops more strongly with T c ff than given by relation (||). 
In case of an amplification, rapid changes of the stellar surface 



parameters would occur, leading to regions in the HRD devoid 
of stars. 

Clearly, we do not observe such "gaps" in the region of the 
HRD studied here, and indeed one can verify with the aid of 
Fig. U that s env as controlled by the surface conditions exhibits 
a stabilizing behaviour since 



d So 



d Tff 



> 0. 



(3) 



L— const 



Similar relations hold for the atmospheric entropy minimum 
(not shown) and the entropy jump (see Fig. Q). Hence, all 
factors influencing s cnv suggest that the convective instabil- 
ity described above will not be encountered. However, at least 
in principle, it cannot be excluded that under special circum- 
stances a convectively driven runaway might occur. 

5.4. The gravity-darkening exponent 

On the surface of a slowly rotating star the relation between 
local T c ff and gravity g can be approximated by 



Teff oc gf. 



(4) 



Sometimes also referred to as giavity-brighteningi 
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Fig. 5. «mlt for standard mixing-length theory (Bohm-Vitense 1958) with A = qmlt Hp (A: mixing-length, Hp: local pressure 
scale height). The presentation of the data is analogous to Fig. [| 



Lucy (1967) argued that the so called gravity-darkening expo 
nent (3 for stars with convective envelopes is given by 

'd log T cff 



dlogg 



(5) 



Figure [| represents lines s onv = const, so the gravity-darkening 
exponent given by our RHD simulations can be readily de- 
duced from this plot. We find an increase of j3 from 0.07 to 0. 10 
when going from the F- to the K-dwarfs, and a slight decrease 
of (3 with decreasing gravity. Basically, this confirms Lucy's re- 
sult (3 — 0.08, with the novelty that our approach eliminates the 
weakest point of his analysis, namely that MLT provides a rea- 
sonable scaling relation s cnv (T e ff, log g) with constant «mlt- 
Current observations of f3 (see e.g. Alencar & Vaz 1997) are 
consistent with our findings, but show a rather large scatter and 
do not allow a critical test of our results. We therefore refrain 
from a further discussion here. 

5.5. Contrasting ours with other approaches 

The long-standing lack of a reliable theory of convection has 
prompted numerous attempts to remedy the situation. In the 
following we comment on the work of two groups and highlight 
the major differences in the involved physics between their's 
and our approach, hoping to clarify the possible reasons for 
deviating results. 



In a series of three papers Lydon et al. (1992, 1993a,b; 
hereafter LFS) presented a formulation of convective trans- 
port based on results of numerical experiments by Chan & 
Sofia (1989). The idealized numerical experiments were set up 
to describe the generic properties of almost adiabatic convec- 
tion; radiative energy transport was included only as a diffu- 
sive flux computed with constant conductivity. Chan & Sofia 
extracted local statistical relations from their numerical data 
which were subsequently used by LFS to derive an expression 
for the convective energy flux suitable for calculating stellar 
structure models. LFS modeled the Sun as well as the A- and 
B-component of the a Centauri system. They found essentially 
the same «mlt in all cases when they translated their formula- 
tion into an effective mixing-length parameter. This is consis- 
tent with our results, but in our opinion not a strong statement 
since the rather small differences in the surface parameters and 
chemical composition among the three stars makes it difficult 
to detect changes of «mlt- Moreover, there is a principle prob- 
lem with the LFS convection formulation since it relies on nu- 
merical experiments for adiabatic convection. In a follow-up 
project to the work of Chan & Sofia, Kim et al. (1995) found 
significantly different statistical relations among the fluctuating 
quantities in the superadiabatic regime, probably implying that 
larger uncertainties have to be attributed to the original results 
of LFS. 
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Fig. 6. ocmt for the Canuto & Mazzitelli convection theory with A = z+ocmt Hp top (z: distance to the upper (Schwarzschild) 
boundary of the convective zone, Hp yt0 p' pressure scale height at the upper boundary). The acMT -calibration reproduces the 
same underlying s env -distribution as the aMLT-calibration shown in Fig. [|. The presentation of the data is analogous to Fig. ^| 



In CMT the main improvement with respect to MLT was 
the inclusion of a larger spectrum of eddies supposed to ex- 
ist in a turbulent medium. Radiation was only included in a 
MLT-like fashion, and effects due to the compressibility of the 
medium were neglected. In our investigation we have regarded 
the characteristic length scale in CMT as a free parameter, ad- 
justing it to fit the results of our simulations. If we compare the 
calibration of CMT shown in Fig. ^ with standard MLT shown 
in Fig. ||, we find a more extended plateau in the solar vicin- 
ity. The general tendency of CMT to make efficient convection 
more efficient and inefficient convection less efficient results in 
a smoother behaviour relative to MLT. But towards high and 
low effective temperatures, a clear over-compensation occurs, 
leading to a change of qcmt- In MLT as well as in CMT, a sig- 
nificant variation of the free parameter is needed, and no sim- 
ple scaling is apparent. We do not find a qualitative difference 
in their ability to reproduce s cnv from the RHD simulations. 

LFS as well as CMT essentially relate the convective flux 
to the local conditions in the flow. Such an approach is not 
well motivated from the simulation results where convection 
appears to be an extremely non-local phenomenon governed by 
the processes taking place in the thin cooling layer at the stellar 
surface (see Spruit 1997 for a recent discussion) while the bulk 
of the convection zone adapts to these layers. Even though we 
fit our simulation results to local theories, we do not assume 



that convection can be described in local terms since the fit is 
of purely formal nature. Furthermore, both LFS and CMT treat 
the effects of radiative transfer in a rudimentary fashion, de- 
spite the well-known fact that for a quantitative description of 
stellar convection a proper treatment of the radiative transfer is 
crucial. Indeed, much effort goes into the realistic modeling of 
radiative transfer in our RHD simulations. CMT describes con- 
vection within the picture of incompressible turbulence. How- 
ever, it is just the compressibility of the stellar gas that gives 
turbulence in stars its special inhomogeneous character (see 
Nordlund et al. 1997 for a detailed discussion). In contrast, our 
RHD simulations fully include the effects of compressibility. 

6. Concluding remarks 

We presented a calibration of the asymptotic entropy s cnv and 
the corresponding mixing-length parameter qmlt for solar- 
type stars basing on radiation-hydrodynamics models. Despite 
the fact that s onv is all that is needed to construct a stellar struc- 
ture model, it was helpful to translate s cn v into a corresponding 
«mlt since ccmlt proved to be less sensitive to the physical 
and numerical input to the models than s env itself. We gave 
only a description of the numerical results without providing 
an explanation in physical terms — not even on a qualitative 
level — of the scaling behaviour of «mlt which we observe. 
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We shall come back to this interesting issue in the next paper of 
this series where we can look at it from a broader perspective 
by including models of sub-solar metallicity. 

Figure |^ shows a remarkably simple structure if one con- 
siders the complex interplay of fluid flow and radiation which 
governs the dependence of s onv on the stellar parameters. Cer- 
tainly, the simplicity of this dependence is a major reason for 
the relative success of MLT to predict s onv , as evidenced by 
the moderate variation of «mlt in Fig. [|. Looking at relative 
changes ignores the more fundamental problem of fixing the 
absolute value of «mlt, which in practice is done by taking 
recourse to empirical constraints. Our RHD models provide a 
determination of the zero point from first principles. Since in 
MLT important pieces of the physics of convection are miss- 
ing — at least within the present physical interpretation of the 
MLT formulae — our work should not be considered as vali- 
dation of MLT, even though MLT is capable of matching some 
of our simulation results quite well. 

Work is under way to check and apply the hydrodynam- 
ical convection models beyond the comparison with the Sun. 
There are classical procedures which allow the determination 
of the convective efficiency at various locations in the HRD 
(position of the red giant branch, shape of the main sequence, 
evolution of binary stars). A lot of work has already been ded- 
icated to empirical determinations, but conclusions are some- 
times conflicting and no clear picture has emerged yet. We sus- 
pect that systematic uncertainties are actually often larger than 
estimated. E.g. Castellani et al. (1999) emphasize discrepan- 
cies in fitting the main sequence of open clusters which are 
related to the temperature-color transformation and the uncer- 
tainties in «mlt- Clearly, with an independent calibration of 
omlt at hand one can disentangle both effects. For the case 
of globular clusters Freytag & Salaris (1999) have studied the 
effects which are expected from our calibration on the shape of 
the turn-off and the position of the red giant branch. By using 
our calibration the effective temperatures of their evolutionary 
models become essentially unaffected by the uncertainties in- 
herent to MLT. The uncertainties related to the temperature- 
color transformation remain present but one can at least judge 
the internal accuracy of the transformation. The precise HIP- 
PARCOS data — in particular for some open clusters — might 
allow a detailed investigation of these issues. Moreover, he- 
lioseismology proved to be an invaluable tool in the case of 
the Sun, and asteroseismic measurements of the internal stellar 
structure appear to be a promising way to gain further insight. 

Already in 2D, the construction of hydrodynamical model 
grids is computationally demanding. Nevertheless, the step to- 
wards 3D models is desirable, and first calibrations of ccmlt 
from 3D models are now available for smaller stellar samples: 
Trampedach et al. (1997) presented a calibration based on 6 
models in the solar vicinity, Abbett et al. (1997) presented a 
calibration for the Sun. The latest version of the calibration of 
Trampedach et al. (Trampedach 1998, priv. comm.) is consis- 
tent with our calibration except for a systematic offset which is 
caused by differences in the employed low-temperature opac- 
ities; the scaling behaviour is similar as far as it is possible to 



judge from the small set of models. Abbett et al. do not provide 
c*mlt with sufficient precision to allow a critical comparison 
with our results. Own work is in progress to quantify system- 
atic differences between 2D and 3D models. As already men- 
tioned, first results for the Sun give a change of «mlt by +0.07 
with respect to the 2D models. We expect systematic changes 
of similar magnitude for other stars. Hence, it appears unlikely 
that conclusions of this paper will have to be altered substan- 
tially when more 3D results become available. 

Last but not least we want to reiterate the warning that 
the presented calibration of ckmlt is only intended to repro- 
duce s onv and the entropy jump. The detailed temperature pro- 
file of the superadiabatic layers is not necessarily represented 
adequately by an MLT model with our calibrated «mlt (see 
Fig. ||). Moreover, the calibration is not suitable for provid- 
ing the optimum mixing-length parameter for convective stellar 
atmospheres. Preliminary results for the Sun (Steffen & Lud- 
wig 1999) show that matching the emergent radiation field of 
multidimensional simulations with ID standard atmospheric 
models requires a quite different mixing-length parameter than 
is needed for matching the entropy jump. We do not consider 
this as a contradiction; it just indicates shortcomings of MLT 
and the ID idealization. 
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Appendix A: Mixing-length formulations 

Generically, the well-known equations of MLT (cf. Cox & 
Giuli 1968) for the convective efficiency 



V - V c 



pc p v c 



V c - V ad / 3 aT3 
the convective velocity 



i+4 

T 



f: 



AV(V - V ) 



Hp 

and the convective flux 

pc p v c TA{\7 - V c ) 



F c = h 



Hp 



(A.l) 



(A.2) 



(A.3) 



contain four dimensionless free parameters fx, f%, fs, and f/±. 
In the equations, the optical thickness of a convective eddy t 
and the isobaric expansion coefficient 5 are defined as 



XPA, 6 



dlnp 
dlnT 



(A.4) 



Further quantities entering the equations are the mixing- 
length A, specific heat at constant pressure c p , Stefan- 
Boltzmann's constant a, density p, temperature T, actual 
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Table A.l. Constants /; for the MLT formulations of various 
authors. (ML1 and ML2 are commonly used terms in the white 
dwarf community.) 



Table B.l. Thermodynamic quantities as a function of entropy 
and pressure from the RHD EOS for solar metallicity and Y = 
0.28. Unless noted otherwise the values are given in cgs-units. 
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temperature gradient V, temperature gradient of a convec- 
tive eddy V c , and adiabatic gradient V a d (V stands for 
d\nT/dlnP). Table Al provides the f-parameters for vari- 
ous MLT formulations. The values of «mlt given in this pa- 
per refer to the MLT formulation by Bohm-Vitense (1958) as 
specified in the first row of the table. 

A fine point in our implementation of MLT concerns the 
evaluation of the radiative gradient V la d- We do not use the 
commonly adopted expression for the radiative gradient in dif- 
fusion approximation, but instead derive V ra( j by differenti- 
ation of the T(r)-relation. As a starting point, consider the 
T(r)-relation for a grey atmosphere in radiative equilibrium 



T 4 = -T c %[r + q(r)} 



(A.5) 



While in the case of a grey atmosphere q is the Hopf-function, 
we can introduce a modified q such that Eq. A.5 represents any 
prescribed T(r)-relation. A reasonable q should have the prop- 
erty 
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8.0 


4.531 


-4.616 


1.415 


1.320 


6.370 


2.00 


9.0 


4.859 


-3.966 


1.534 


1.081 


4.169 


2.00 


10.0 


5.201 


-3.325 


1.646 


1.026 


3.520 


2.00 


11.0 


5.598 


-2.723 


1.666 


1.004 


3.352 


9 no 


1 9 O 




-Z. LAJ 


1 .oo / 


1 nrn 




2.00 


13.0 


6.398 


-1.523 


1.667 


1.002 


3.342 


2.25 


8.0 


4.789 


-4.899 


1.474 


1.090 


4.641 


2.25 


9.0 


5.124 


-4.249 


1.655 


1.010 


3.419 


2.25 


10.0 


5.523 


-3.648 


1.667 


1.001 


3.341 


2.25 


11.0 


5.923 


-3.048 


1.667 


1.000 


3.339 


2.25 


12.0 


6.323 


-2.448 


1.667 


1.000 


3.338 


2.25 


13.0 


6.723 


-1.848 


1.667 


1.000 


3.338 


2.50 


8.0 


5.049 


-5.174 


1.660 


1.004 


3.377 


2.50 


9.0 


5.448 


-4.574 


1.667 


1.000 


3.338 


2.50 


10.0 


5.848 


-3.974 


1.667 


1.000 


3.338 


2.50 


11.0 


6.248 


-3.374 


1.667 


1.000 


3.338 


2.50 


12.0 


6.648 


-2.774 


1.667 


1.000 


3.338 


2.50 


13.0 


7.048 


-2.174 


1.667 


1.000 


3.338 



which guarantees that V ra d becomes independent of r and con- 
sistent with the diffusion approximation for large optical depth 
(provided the Rosseland scale is adopted). The advantage of 
this procedure is that we obtain a smooth and realistic V ra d 
which is consistent with the specified T(r)-relation irrespective 
of optical depth. We never have to distinguish between atmo- 
sphere and optically thick layers and can use MLT throughout, 
for the small price of having to compute the r-scale as well. 

Appendix B: Entropy scale 

The entropy is not always available from an EOS, and if so it 
involves an arbitrary additive constant. In order to relate en- 
tropy values given in this paper to more common thermody- 
namic variables, we provide a small list of quantities as a func- 
tion of entropy and pressure in Table Bl, calculated from the 
EOS which was used in the RHD simulations. We have con- 
centrated on higher temperatures characteristic of the deeper 
layers of solar- type stellar envelopes. The data should allow to 
interrelate values for s cn v to pressures, temperatures, and den- 
sities, and provide the entropy zero point of our EOS. For in- 



terpolation purposes one might take advantage of the following 
differential thermodynamic relations 



1 ds 
d In p = — d In P — 8 — . 

Ti c p ' 



and 



d\nT = V ad dln.P 



ds 



(B.l) 



(B.2) 



V a d is related to the thermodynamic derivatives given in Ta- 
ble B 1 according to 



V ad = 



PS 
pc p T' 



(B.3) 



Appendix C: Fitting functions 

The fits for s cnv , entropy jump As, q;mlt, and «cmt s how n 
in Figs. §, @J|, and | were computed from expressions dC3j), 
(C.4), and (jC.5) (for omlt an d ocmt) in terms of the auxiliary 
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Table C.l. Coefficients for fitting functions 



coefficient 


Scnv / SO 


As/so 


o;mlt 


QCMT 


d() 


1.6488 


0.029 


1.587 


0.414 


ai 


0.0740 


0.166 


-0.054 


-0.100 


0,2 






0.045 


-0.026 


<13 


1.7860 


1.225 


-0.039 


-0.065 


(24 


-1.6762 


-1.209 


0.176 


-0.092 


as 


0.1274 


0.073 


-0.067 


0.193 


a@ 


-0.1412 


-0.132 


0.107 


-0.125 



variables T and g denned below. The coefficients Oj are listed 
in Table CI. The fitting functions were generated automatically 
from a prescribed set of basis functions with the mere intention 
of providing numerical fits to the RHD data points which then 
allow the RHD results to be utilized easily in other contexts. 
Physical considerations did not enter, and we present the for- 
mulae in a form more suitable for computing rather than for 
human interpretation. We warn the reader that the fits quickly 
loose their meaning outside the region covered by the RHD 
model grid and should not be extrapolated too far. 

f = (T cS - 5770)/1000 (C.l) 

g = log(.g/27500) (C.2) 

= a + a 5 f + a e g + a\ exp[a 3 T + a 4 g] (C.3) 

So 

As ~ 

= a + Oi cxp[(a 3 + a 5 T + a 6 g)T + a 4 g] (C.4) 

so 

a = a + (a! + (a 3 + a 5 f + a 6 g)f + a 4 g)f + a 2 g (C.5) 
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